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Is Newton's gravity sufficient to handle the weakly nonlinear evolution stages of the cosmic large- 
scale structures? Here we resolve the issue by analytically deriving the density and velocity power 
spectra to the second order in the context of Einstein's gravity. The recently found pure general 
relativistic corrections appearing in the third-order perturbation contribute to power spectra to the 
second order. In this work the complete density and velocity power spectra to the second order are 
derived. The power transfers among different scales in the density power spectrum are estimated 
in the context of Einstein's gravity. The relativistic corrections in the density power spectrum are 
estimated to be smaller than the Newtonian one to the second order, but these could be larger than 
, higher-order nonlinear Newtonian terms. 

■ PACS numbers: PACS numbers: 98.80.-k, 04.25.Nx, 98.80.Jk 



oo 

o 

O 



6 ■ 

H ■ 



I. INTRODUCTION 



The weakly nonlinear process of gravitating system has fundamental importance in cosmology [U, 0) S 0|- The 
large-scale structure in the observable universe apparently show its nonlinear nature. Considering the success of the 
Friedmann cosmology with its spatial homogeneity and isotropy assumptions, it is likely that in near horizon scale 
the structures are in near linear stage. In the small scale, however, the cosmic structure (say, distribution and motion 
of galaxies) are in fully nonlinear stage. In between the two scales, we have weakly nonlinear stage. Even for the 
nonlinear structures in the present epoch, as the cosmological structure grows under gravity from linear to nonlinear 
stages there must be transition era which can be regarded as weakly nonlinear. The density and velocity power 
: spectra provide the main observations of the large-scale structure which can be directly compared with theories of the 
structure formation. All theoretical studies of the weakly nonlinear evolution of the large-scale structure have been 
based on Newton's gravity [U, S H]- Newton's gravity sufficient to handle the situation? In this work we will 
' resolve the issue by analytically deriving the power spectra to the second order in the context of Einstein's gravity. 
00 We will derive the second-order density and velocity power spectra in the context of relativistic cosmology for the 

first time. In a zero-pressure medium the pure general relativistic corrections appear in the third order; this was 
^ ' only recently shown by us in The third-order relativistic corrections, however, contribute to the density and 
• i-H , velocity power spectra even to the second order which was unknown previously. We show that compared with the 
Newtonian contributions to the second-order power spectra the relativistic contributions are generically multiplied by 
^ . a factor (scale/horizon scale) squared, thus suppressed in the small scale. In the Newtonian theory, in the second- 
■ - - ' order density power spectrum, the nonlinear power transfer from large-scale nonlinearity is known to exactly cancel 
to the leading order, thus opening a possibility of relativistc effect becoming important. In this paper we will show 
that, even in such a situation, the pure general relativistic contribution to the density power spectrum is smaller 
than the remaining second-order Newtonian ones. We conclude that, even in the context of Einstein's gravity, up 
to the second order, the fc* long wavelength tail in the density power spectrum previously known in the Newtonian 
analysis is the only important effect of the nonlinear power transfer among different scales. However, the leading 
order relativistic contribution in the second-order power spectra could be larger than the higher (third and higher) 
order pure Newtonian contributions to the power spectra , thus demanding a caution in the pure Newtonian study. 
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II. BASIC EQUATIONS 

We consider an irrotational dust (zero-pressure fluid) without the gravitational waves in a flat Friedmann back- 
ground. To the third order the basic perturbation equations in Einstein's gravity are recently derived in [5]. These 
are 
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where d and u are the density perturbation [S = Sg/g) and the perturbed velocity, respectively; a is the cosmic scale 
factor and g is the background density. A perturbed order variable 1^9 is a metric perturbation variable related to the 
perturbed three-space curvature. To the linear order (p can be related to 6 and u as Q 



[S^ + aA-^V ■ u) , 



(4) 



where (5$ is the perturbed order Newtonian gravitational potential; 6^ is related to S by Poisson's equation 

A$ = AnGga^S. 



(5) 



Up to the second order, remarkably, Eqs. ([T]) and Q coincide exactly with the ones in Newtonian theory; we call 
it a relativistic/Newtonian correspondence to the second order Q. This is why we simply call 6 and u as the 
(Newtonian) density and velocity perturbations even in the present relativistic situation; in the relativistic context S, 
u and if are related to the certain gauge-invariant combinations of variables, see [1, l^. It is also remarkable to notice 
that all the pure relativistic third-order correction terms have ip factor compared with the relativistic/Newtonian 
second-order terms. Contribution of the gravitational waves to the third order can be found in [5]; for equations 
in the multi-component case, see the third reference in 0]. We note that as the above equations are derived in the 
relativistic perturbation theory, these are valid in the fully general relativistic situation and in all scales (including 
the super-horizon scale) as long as the perturbative assumption is met. 



III. SOLUTIONS IN THE PHASE SPACE 



As we are considering a flat background, we may take the Fourier transformation deflned as F(k) — J d'^xF{x.)e^^'^. 
By introducing u = Vu and = V • u = Aw, we have u(k, t) = — iku(k, t) and 0(k, t) = —k'^u{'k,t). In the phase 
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space, Eqs. (H])-® become 
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where C+ indicates terms replacing k' to k — k'; C++ indicates two sets of terms, one replacing k' to k — k' — k", 
and the other replacing k" to k — k' — k". 

In order to derive the perturbative solutions we expand 



5{k,t) ^ 5i{k,t) + 52{k,t) + 5^{k,t) 



(9) 



and similarly expand 8 and Lp. To the linear order, Eqs. (O and ([7]) give 5i + 2{d/a)5i — AirGgSi = 0. Equations 
up to this point are valid in the presence of the cosmological constant. In the following, in order to derive analytic 
solutions we assume an absence of the cosmological constant. In a flat background without the cosmological constant, 
we have a oc t'^^^, GnGg ~ t^"^ , thus we have two solutions i5i(k, t) oc t^/^ and We ignore the decaying solution in 
an expanding phase and set 

5i{k,t) = A{k)e"^^^h'^/^, Bi{k,t)^--A{k)e"f'^^^at-^'^, ipi{k,t) ^- {t/lnf t^''^A{k)e"l'^^\ (10) 
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In a flat background without the cosmological constant, to the linear order, from Eqs. ^ and ([5]) we have ip = 
— (5/3)5$/c2. We introduced a scale £ = a/k and the Hubble horizon scale £h = c/{a/a); thus, = a/{kc). 

Notice that ipi is time independent. To the second order from Eqs. ([6]) and ([7]) we have the solutions 
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J(k, k', k - k') s [4F(k, k', k - k') + hH{k, k') + 5H{k, k - k')] , 
L(k, k', k - k') = [8F(k, k', k - k') + 3i7(k, k') + 3i7(k, k - k')] , 
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In order to derive third-order solutions we need X to the second order. Using Eq. ^TU\i . Eq. ([5]) becomes 
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To the third order from Eqs. ([6]) and using Eqs. ((T0l) - ((T5)) . we have the sohitions 

AF(k, k', k - k')L(k - k', k", k - k' - k") 

\H{k, k') J(k - k', k", k - k' - k") + \H(k, k - k')L(k - k', k", k - k' - k") 
18 18 



t^^±_ J Sk' j d3fc''A(k'0e''^('''')|A(k0A(k-k'-k'0e'^(''')+*"*(''-'''-'''' 

e / (k-k^-k")-k" 4 (k' + k") • k" (k - k' - k") ■ k" k ■ (k - k^ - k'^) 
fc^^2 ^3 fc^^2 |k-k' -k"|2 



A(k')A(k - k' - k")e*'^(''')+*'^(''-'''-''"' Af (k - k', k", k - k' - k")' 



fc2 



k-k' 



/I 3k -k' 15 k' -(k-k') 3 fc2 k' • (k-k')\ + 
^ 1^2 ~ 21^2" |k-k'|2 +2fc^ |k-k'|2 j + 

6l3(k,i) = -^at^^ y" d^fc'y" d3fc"A(k")e"^(*'")|A(k')A(k-k'-k")e*"*(''')+*'^(''"'''-''" 
F(k, k', k - k')i(k - k', k", k-k' - k") 



-iiJ(k, k') J(k - k', k", k-k' - k") + iiJ(k, k - k')L(k - k', k", k - k' - k 



/ '''7 ^'fc"^(k")e^^^'''7^(k'M(k-k'-k") 



i0(k')+i0(k-k'-k") 



fc2 / 2 (k - k' - k") ■ k" 16 (k' + k") • k" ^ 16 (k - k' - k") ■ k" k • (k - k' - k"; 
^ l~3 ^ 



fc"2 



9 



fc"2 



A:"2 



|k-k'-k"|2 



(16) 



C++ 



A(k')A(k - k' - k")e*"^"'')+*'^(''-'''-''")M(k - k', k", k-k' - k") 



fc2 



|k-k' 



2 k-k' 3 k' • (k - k') fc2 k' • (k - k') 

^ ' "3^^~fc^"2 |k-k'|2 "^A^ |k-k'|2 



(17) 



We note that the pure general relativistic contributions first appearing in the third order are generally multiplied by 
a [J-IIh)'^ factor which came from terms in Eqs. (H])-©, see Lp in Eq. ([10]); {IjlH)^ is small in the small scale but 
becomes of order unity in near horizon scale. The Newtonian part of ^3 is proportional to a?" , and the relativistic 
part is proportional to . The Newtonian part of ^3 is proportional to ai, and the relativistic part is proportional to 



IV. POWER SPECTRA 



The density power spectrum is 
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where indicates the real part. The velocity power spectrum |^(k, t)|2 can be similarly expanded. Assuming the 
random phase, the second term in the right-hand-side vanishes. The second-order power spectra of density and 
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velocity follow from Eqs. (fT0l) - ((T7)) . These are 
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Notice that, in the second-order power spectrum, compared with the Newtonian contributions the pure general 
relativistic effects are simply multiplied by a {£/£h)'^ factor. Newtonian part of the density spectrum is proportional 
to a* and the relativistic part is proportional to a^. In the case of the velocity power spectrum, Newtonian part is 
proportional to and the relativistic part is proportional to . Newtonian parts of the power spectra are known in 
the literature 0, |3| - The pure relativistic contributions to the second-order power spectra are our new contribution 
in this work. 

(i) For the power transfer from the small-scale, thus k' — > oo, we have (in the following we assume isotropic power 
spectrum, thus |(5(k, t)| = |(5(A;,i)|, etc) 
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Using \S{k,t)\'^k^ , with Sk a density contrast at a given wavenumber fc, and introducing a small-scale cut-off 

frequency kc we have 



\Sik,t)\' ^Mk,t)\'+c,{k/k,)^k-'6t+c2Mk,t)\'{k/k,yst+C3Mk,t)\'i£/£HrSi 
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where q are constants of order unity. The ci term is the well known fc"* long wavelength tail generated by the power 
transfer from the small-scale nonlinearity; this was first shown by Zel'dovich [l„i2,, jjj. The Newtonian and relativistic 
contributions in 2TZ{51S3) are smaller than the linear term by factors {k/kcYSj. and (£/^^r)2(5| , respectively; we have 
k/kc = £c/£ < 1 with 4 = a/kc. 

(ii) For the power transfer from the large-scale, thus k' 0, we have 
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2l7r2 



k^\5i{k,t)f / dk'\Si{k',t)\^ 



(23) 



6 



To the lowest order in k' 0, the Newtonian contributions, which have order fc^|(5i(fc, t)p J dk'\Si{k' ,t)\'^ , exactly 
cancel out to the second order; this was shown by Vishniac [1, The next order terms in that limit have order 
\5i{k,t)\'^ J k''^dk'\6i{k' jt)]"^ . Thus, using a large-scale cut-off frequency kc we have 

\S{k,t)\^ ^Mk,t)\' + c,\SAk,t)\Hl + c5\6Ak,t)\'{£/e^^^^ (24) 

Thus, the second-order terms, both Newtonian and relativistic, are smaller than the linear ones by factors S'^ 
and {t/ tuYik/kcYSl. , respectively. Compared with the Newtonian second-order term the relativistic one has 
~ (^/^ff )^(fc/fcc)^ ^ {^c/^hY which is smaller than unity but has no scale dependence. 

V. DISCUSSION 

Our results show that even to the second order the power spectra receive corrections from pure relativistic third- 
order perturbations. Compared with the Newtonian terms, the relativistic contributions are multiplied by a factor 
{(-/ (■hY' , thus generally suppressed in the small-scale limit, but comparable in near horizon scale. In near horizon scale, 
however, as the perturbations are supposed to be in near linear stage, the second-order contributions arc negligible 
compared with the linear-order terms. In the case of the power transfer from the large-scale nonlinearity it was 
previously known that the leading order Newtonian nonlinear contribution in the density power spectrum cancels 
exactly. Even in that situation our investigation shows that the relativistic effect is subdominant compared with the 
remaining Newtonian nonlinear effect. Our analysis of the asymptotic cases shows that the fc^ long wavelength tail 
of density power spectrum generated by the small-scale nonlinearity in the Newtonian theory is the only important 
effect of the nonlinear power transfer even in Einstein's gravity. 

Resolution of the issue of whether the relativistic contributions are always smaller than the Newtonian nonlinear 
effects requires quantitative estimation of the general power spectra presented in Eqs. (|19p and (j20p . This may depend 
on the specific form of linear density power spectrum, and may require numerical integration of Eqs. ^9]) and ([20)l . 
As the higher perturbational order contribution to the power spectrum is in general suppressed by power of (5^ term 
it is likely that such Newtonian contributions are smaller than the relativistic contribution to the second order; 
quantitative estimation may depend on the linear spectrum and the scale. Although the Newtonian perturbation 
theory has a recursion formulae to all orders Q , the relativistic situation should be handled at each order separately. 
Estimation of non-Gaussian contribution of the pure relativistic corrections is also an important and interesting issue 
to be addressed: our Eqs. (I6|)-([8|) or Eqs. (|10 p - (fT7|l provide the starting point for evaluating the bispectrum or the 
higher order correlation functions. In the presence of the cosmological constant we have to go back to Eqs. ([HI)-® 
which are valid in the presence of the cosmological constant; for a general expression of Lp in terms of 5 and 6 to the 
linear order, see Eq. ([¥]). Our perturbation equations are fully relativistic while assuming perturbations to be weakly 
nonlinear. In the small scale where structures are in fully nonlinear stage while the relativistic effects are small, the 
cosmological post- Newtonian approximation provides a complementary theoretical framework to handle the structure 
formation f^. Density power spectrum based on the cosmological post-Newtonian equations has not been studied in 
the hterature. Investigation of these issues are left for future studies. 
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